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REMARKS/ARGUMENTS 

With this amendment, claims 53, 55-58 and 60-74 are pending. For convenience, 

the Examiner's rejections are addressed in the order presented in an October 28, 2005 Office 
Action. 

L Rejections under 35 U.S.C, §112, first paragraph, written description 

Claims 53, 55-58 and 60-74 are rejected under 35 U.S.C. §112, first paragraph for 
allegedly failing to comply with the written description requirement. According to the Office 
Action, the specification lacks description of the claimed invention, such that a skilled artisan 
would recognize that Applicants had possession of the claimed invention at the time of filing. In 
particular, the Office Action appears to assert that the genus of glycosyltransferases, accessory 
enzymes, and oligosaccharides are not adequately described because amino acid sequences of 
known proteins are not included in the specification. Applicants respectfiilly traverse the 
rejection. 

The claims are directed to cell-based methods of producing product saccharides 
by combining a heterologous accessory enzyme for forming a nucleotide sugar and a 
heterologous glycosyltransferase in the same cell; contacting the cell with an acceptor saccharide 
and allowing formation of the nucleotide sugar through activities including the activity of a 
heterologous accessory enzyme, followed by transfer of the sugar to the acceptor saccharide 
through the activity of the heterologous glycosyltransferase. 

The Office Action cites Regents of the University of California v. Eli Lilly, 43 
USPC2d 1398 (Fed. Cir. 1997) as standing for the proposition that, in order to describe a genus, 
AppHcants must provide at least one species of the claimed genus and identify the common 
characteristics of the claimed molecules. However, a recent opinion fi-om the Federal Circuit of 
Appeals reviewed the description requirement for known sequences and has distinguished Lilly, 
See, e,g„ Capon v. Eshar, 76 USPQ2d 1078 (Fed. Cir. 2005). 

The dispute in Capon began with an interference proceeding; and the claims were 
directed to a chimeric gene that encodes a fiision between an endogenous protein expressed on 
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the surface of immune cells and a variable region of a antibody. Neither party provided a full 

length DNA or encoded amino acid sequence of a chimeric gene. The Board of Patent Appeals 

and Interferences ruled on the patentabiHty of claims presented by both parties and, citing ruling 

precedent, e.g., Lilly, ruled that both claim sets were invalid because the chimeric genes were 

described only in terms of function and a sequence of a chimeric gene was not included in the 

application. Capon, 76 USPQ2d at 1082. 

The Federal Circuit found that the rule put forth by the Board was "an 

inappropriate generalization", stating: "When the prior art includes the nucleotide information, 

precedent does not set a per se rule that the information must be determined afresh." Capon, 76 

USPQ2d at 1084-1085. The Federal Circuit made a distinction over Lilly, because, unlike the 

claims of Lilly, the claims at issue in Capon did not concern a discovery of gene function or 

structure and thus, sequence information did not add description. Capon, 76 USPQ2d at 1085. 

The chimeric genes here at issue are prepared from known DNA 
sequences of known function. The Board's requirement that these 
sequences must be analyzed and reported in the specification does 
not add descriptive substance. The Board erred in holding that the 
specifications do not meet the written description requirement 
because they do not reiterate the structure or formula or chemical 
name for the nucleotide sequences of the claimed chimeric genes. 

As in Capon, the present claims are directed to methods to use known proteins to 
make known compounds. The present claims are not directed to the discovery of new protein or 
oligosaccharide functions or structures. Accessory enzymes, glycosyltransferases, and product 
saccharides are adequately described in the specification and the description enables those of 
skill in the art to use the claimed methods. Moreover, the specification provides representative 
numbers of each genus, i.e., accessory enzymes, glycosyltransferases, and product saccharides. 
For example, the specification provides specific examples of at least 12 accessory enzymes for 
forming a nucleotide sugar. Specification at pages 29-34. Table 1 at page 28, also provides nine 
examples of cycles for forming nucleotide sugars, including the enzymes required to regenerate a 
nucleotide triphosphate from a nucleotide. Using a single enzyme as an example, at page 32 of 
the specification, ten UDP-Gal pyrophosphorylase enzymes with accession numbers are 
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disclosed. Six are from prokaryotic or viral organisms: Lactobacillus^ E. coli, B. subtilis. 
Neisseria, Haemophilus, and Streptomyces. Four are from eukaryotes: Rat, bakers yeast, mouse, 
and human. Similar analysis is done for the other enzymes disclosed throughout 6 pages of the 
specification. Thus, by providing detailed information about the identity of multiple accessory 
enzymes in the specification combined with well-known characteristics of enzymatic pathways 
to synthesize nucleotide sugars, applicants have adequately described of this genus. 

Glycosyltransferases are similarly well described in the specification. Eukaryotic 
glycosyltransferases are described from page 19, line 1 through page 21, line 15. Embodiments 
include fixcosyltransferases, galactosyltransferases, sialyltransferases, glucosyltransferases, N- 
acetylgalactosaminyltransferases, and mannosyltransferases. References are cited for enzymes 
from many different organisms and some Accession Numbers are given. For example, 
galactosyltransferases are listed starting at page 19, line 25 through page 20, linel4. Twelve 
different galactosyltransferases are disclosed from human, bovine, murine, and porcine sources, 
as well as an enzyme form the yeast Schizosaccharomyces pombe, a(l,3), a(l,2), and q(1,4) 
galactosyltransferases are disclosed, as is a ceramide galactosyltransferase. Similar disclosure is 
made for the other eukaryotic glycosyltransferases. 

Many prokaryotic glycosyltransferases are also disclosed starting at page 21, line 
16 through page 22, line 25. Prokaryotic glycosyltransferases include enzymes involved in 
synthesis of LOS. Disclosed prokaryotic enzymes include galactosyltransferases, 
glucosyltransferases, N-acetylglucosaminyltransferases, fijcosyltransferasese, and 
glycosyltranseferases involved in the synthesis of structures containing lacto-N-neotetraose. The 
enzymes come from a variety of prokaryotic sources, including E. coli. Salmonella typhimurium, 
Klebsiella pneumoniae, S, enterica. Yersinia enterocolitica, Mycobacterium leprosum, 
Pseudomonas aeruginosa. Neisseria gonorrhoeae, and A^. meningitidis. References are cited for 
each enzyme disclosed and some accession numbers are given. Thus, the description of 
glycosyltransferases used in the claimed methods is sufficient to meet the written description. 

Both accessory enzymes and glycosyltransferases are well known to those of skill 
in the art. In addition to the disclosure found in the specification, numerous reference materials 
are available that discuss accessory enzymes and glycosyltransferases, for example biochemical 
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textbooks. See e.g. Lehninger, Principles of Biochemistry (1984); Stryer, Biochemistry (1995); 
Zubay, Biochemistry (1986). As discussed above, neither accessory enzyme, nor 
glycosyltransferases must be described in detail in the specification. Hybritech, 802 F.2d at 1384. 

The Office action also appears to assert that the description of oUgosaccharide 
products in the specification is not adequate. Applicants respectfixUy point out that, 
oligosaccharide products were well known at the time of filing. Known compounds are not 
required to be described in detail. The application provides adequate description of 
oligosaccharides, including structural information. For example, product oligosaccharides are 
described at Table 1, page 45; at page 48, Unes 19-25; at page 49, line 1 through page 5, line 10 
(including Table 2); at page 52, Hne 10 through page 53, line 28; and at page 54, line 25 through 
page 57, line 10 (including Table 3). Thus, the description of product oligosaccharides made 
using the claimed methods is sufficient to meet the written description requirement. 

Based on the disclosure of the specification, one of skill would be able to practice 
the claimed methods. The Office Action's insistence that the sequence of accessory enzymes and 
glycosyltransferases must be disclosed in the specification to meet the description requirement is 
not appropriate under current Federal Circuit case law. 

The Office Action also responded to AppUcants' discussion in the previous 
response of Example 18 of the USPTO's Synopsis of Application of Written Description 
Guidelines. According to the Office Action, the subject matter of the pending claims differs 
fi-om that of Example 18 because the pending claims encompass genera with allegedly 
"substantial variation between members." Applicants respectfully remind the Examiner that the 
claims of Example 18, like the pending claims, base their novelty on the steps of the claimed 
methods. The language of the adequately described claim is as follows: "A method of 
producing a protein of interest comprising; obtaining Neurospora crassa mitochondria. . . .and 
recovering said protein of interest." Emphasis added. Thus, the adequately described protein in 
Example 18 is not described in terms of function, not described in terms of structure, but is only 
described in terms of the users subjective interest. Applicants respectfully ask the Examiner to 
specifically explain how a genus of "proteins of interest" does not encompass "substantial 
variation between members." The analysis of the method claimed in Example 18 was performed 
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on the method steps, not the recited "protein of interest." Example 18 did not require listing of 
perhaps millions of amino acid sequences that could be expressed in Neurospora crassa 
mitochondria to meet the written description requirement. Reduction to practice of a single 
embodiment was sufficient. Applicants respectfully request that the Examiner similarly analyze 
the method steps of the pending claims with regard to the written description requirement. 

In view of the above arguments, withdrawal of the rejections for alleged lack of 
written description is respectfully requested. 

IL Rejections under 35 U.S.C. §103(a) 

Various combinations of claims are rejected for alleged obviousness in view of 
various combinations of references. Applicants respectfully traverse the rejections. 

The Office Action has not estabhshed a case of prima facie obviousness. To 
estabUsh a case of prima facie obviousness, the Examiner must meet three basic criteria: 

First, there must be some suggestion or motivation, either in the references 
themselves or in the knowledge generally available to one of ordinary skill in 
the art, to modify the references or to combine reference teachings. Second, 
there must be a reasonable expectation of success. Finally, the prior art 
reference (or references when combined) must teach or suggest all the claim 
limitations. In re Vaeck, 947 F.2d 488, 20 USPQ2d 1438 (Fed. Cir. 1991). 
M.P.E.P. §§ 706.02(j) and 2143. 

The references cited by the Examiner fail to provide a reasonable expectation of 
success in practicing the invention and fail to provide a motivation for the combination of the 
references. In addition, the references cited by the Examiner fail to provide all the elements of 
the rejected claims. Therefore, Applicants respectfully traverse the rejections. 

Claims 53, 56-58, 60, and 65-74 are rejected as allegedly obvious over Samain et 
al, Carbohy. Res. 302:35-42 (1997) in view of Ulbich and van Putten, /. Bacterial. 177:6902- 
6909 (1995). Applicants respectfully traverse the rejection. The claimed invention is a cell 
based method of producing an oligosaccharide, by using a host cell that comprises a 
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heterologous glycosyltransferase and a heterologous enzymatic system to produce a sugar 
nucleotide. The cited references, either alone or combined, fail to provide all the elements of the 
claimed invention and fail to provide motivation for their combination to arrive at the claimed 
invention. 

Samain et al. discloses synthesis of chitoohgosaccharides (COs) in E. coli by 
expression of the A, caulinodans NodC gene, which encodes a CO synthase. Samain et al does 
not disclose or suggest co-expression of an accessory enzyme for enzymatic synthesis of 
GlcNAc, one of the substrates of the NodC protein. In fact Samain et al teach av^ay from use of 
anything other than the endogenous E. coli sugar nucleotide pools, stating at page 36, left column 
that UDP-GlcNAc is maintained at "quite high" intracellular levels in most bacteria, thereby 
providing ample substrate for the NodC protein. Samain et al. also caution that over expression 
of gene products can be unfavorable to oligosaccharide synthesis by causing "unnecessary 
metabolic buming, a growth inhibition by overproduced proteins, and problems of precipitation 
and misfolding of proteins." Samain et al, page 39, left column. Thus, Samain et al wam those 
of skill against over expression of proteins generally, including over expression of multiple 
proteins. 

Ulkich and van Putten discloses cloning of GlcN Ac- 1 -phosphate uridyltransferase 
(GlrnU) gene from gonococcus. Ullrich and van Putten does not disclose or suggest co- 
expression of a glycosyltransferase for enzymatic synthesis of an oligosaccharide. Ullrich and 
van Putten do not discuss or investigate the affect of GlmU expression on levels of either UDP- 
GlcNAc or oligosaccharides containing GlcNAc. Given Samain et al's disclosure that UDP- 
GlcNAc levels are quite high in most cells, one of skill would have no reason to suppose that 
expression of GlmU would have any effect on cellular UDP-GlcNAc levels and therefore would 
not have reason to attempt to co-express, e.g., NodC and GlmU in the same cell. Therefore, the 
combination of Samain et al and Ulhich and van Putten does not render the claims obvious. 

Claims 61-65 are rejected as allegedly obvious over Samain et al in view of 
Ullrich and van Putten, and in ftirther view of Bulow et al Trends Biotech, 9:226-231 (1991). 
Applicants respectfully traverse the rejection. The claimed invention is a cell based method of 
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producing an oligosaccharide, by using a host cell that comprises a fusion of a heterologous 
glycosyltransferase and a heterologous accessory enzyme to produce a sugar nucleotide. The 
cited references, either alone or combined, fail to provide all the elements of the claimed 
invention and fail to provide motivation for their combination to arrive at the claimed invention. 

The above discussion of Samain et al and Ulbich and van Putten is applies to this 
rejection, but will not be repeated here. This analysis discusses the citation of Bulow et al. 
Bulow et al. disclose generally construction of bi-functional proteins by fusing proteins with 
different functions together. Bulow et ai does not discuss or suggest fusion proteins that include 
a glycosyltransferase and an accessory enzyme. 

The Office Action alleges that Bulow et aL discloses favorable enzyme kinetics 
through fusion of two enzymes, /.e., a so-called proximity effect. Bulow et al. was published in 
1991, and by the time the priority application was filed in November 1998, the proximity effect 
had largely been discredited by those of skill. 

First, researchers established that the proximity effect cannot be explained by 
simple diffusion between active sites on the same protein and proposed the existence of 
electrostatic regions on the surface of the fusion protein that allowed substrate channeling 
between the active sites of the enzyme components. See, e.g., Elcock and McCammon, 
Biochemistry 35:12652-12658 (1996); submitted herein as Exhibit A. 

Researchers then began to doubt the substrate channeling proximity effect. 
Naturally-occurring bi-fiinctional fusion proteins that supposedly exhibit the proximity effect had 
been identified and their structure had been solved. See, e.g., Trujillo et al. Prot. Eng. 10:567- 
573 (1997); submitted herein as Exhibit B. Using the structure of a naturally occurring bi- 
fimctional fusion protein as a model, Trujillo et al. constructed a synthetic bi-functional fusion 
protein with mono-functional proteins from a different organism. However, even with kinetic 
and structural information of the naturally occurring bi-functional enzyme available, Trujillo et 
aL were unable to demonstrate a proximity effect in the synthetic bi-functional fusion protein. 
Thus, after Bulow et al. was pubhshed, researchers first found it necessary to revise the 
proximity effect theory to fit new data, and then found that the proximity effect theory could not 
be verified, even under the most optimal conditions, i.e., with structural information supposed to 
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allow construction of a bi-functional protein that exhibits the proximity effect. Thus, at the time 
of filing, based on additional analysis of the proximity effect in e.g., Elcock and McCammon, 
and Trujillo et al., those of skill would not have been motivated to combine the cited references 
to arrive a -bi-functional fusion protein comprising an accessory enzyme and a 
glycosyltransferase, and would not have had an expectation of success in doing so. 

Trujillo et al was pubhshed in May 1997, well before the November 1998 
priority date. Mmore papers have been published discrediting the proximity effect, including a 
paper listing Leif Bulow of Bulow et al as an author. See, e.g., Pettersson et al Eur. J, 
Biochem, 267:5041-5046 (2000), submitted herein as Exhibit C. Thus, those of skill would not 
have found motivation to fuse a glycosyltransferase and an accessory enzyme based on the 
supposed disclosure of "favorable enzyme kinetics" by Bullow et al. 

According to the Office Action Bulow et al, also disclose easy purification of the 
fused enzyme. In fact, Bulow et al.'s discussion of protein purification improvements is limited 
to fusion of an "affinity handle" fused to a protein of interest. Bulow et al. at page 230. The 
affinity handle is chosen because it faciHtates use of ion exchange, hydrophobic, covalent, and 
metal chelate separations. Examples included sugar binding proteins, e.g., beta-galactosidase; 
poly-amino acid sequences; and coenzymes. Neither an accessory enzyme or a 
glycosyltransferase is suggested as an affinity handle to improve protein purification. 

Moreover, the claimed methods are cell based methods. Purification of the 
proteins provides no advantage to the user and will likely slow the process of oligosaccharide 
production, greatly reducing efficiency. Applicants respectfully remind the Examiner that the 
product of the claimed methods is an oligosaccharide, not a protein. Therefore, those of skill 
would not have found motivation to fuse a glycosyltransferase and an accessory enzyme based 
on the supposed disclosure of improved protein purification by Bullow et al. 

Bulow et al. provide no additional disclosure of method steps and provide no 
motivation to use a fusion of a glycosyltransferase and an accessory enzyme in the claimed 
methods. Therefore, the combination of Samain et al., Ullrich and van Putten, and Bulow et aL 
does not render the claims obvious. 
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In view of the evidence submitted with this response and the arguments above. 



AppUcants respectfully request that the rejection for alleged obviousness be withdrawn. 



In view of the foregoing, Applicants believe all claims now pending in this 



Application are in condition for allowance. The issuance of a formal Notice of Allowance at an 
early date is respectfully requested. 



If the Examiner believes a telephone conference would expedite prosecution of 



this application, please telephone the undersigned at 415-576-0200. 



TOWNSEND and TOWNSEND and CREW LLP 

Two Embarcadero Center, Eighth Floor 

San Francisco, Cahfomia 941 1 1-3834 

Tel: 415-576-0200 

Fax: 415-576-0300 

Attachments 

BLK:blk 

60758862 v1 



CONCLUSION 



Respectfully submitted, 
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abstract: Brownian dynamics simulations were performed to investigate a possible role for electrostatic 
channeling in transferring substrate between two of the enzymes of die citric acid cycle. The diffusion 
of oxaloacetate from one of the active sites of malate dehydrogenase (MDH) to the active sites of citrate 
synthase (CS) was simulated in die presence and absence of electrostatic forces using a modeled structure 
for a MDH-CS fusion protein. In the absence of electrostatic forces, fewer than 1% of substrate molecules 
leaving the MDH active site are transferred to CS. When electrostatic forces are present at zero ionic 
strength however, around 45% of substrate molecules are successfully channeled. As expected for an 
electrostatic mechanism of transfer, increasing the ionic strengdi in the simulations reduces the calculated 
transfer efficiency. Even at 150 mM however, the inclusion of electrostatic forces results in an increase 
in transfer efficiency of more than 1 order of magnitude. The simulations therefore provide evidence for 
the iiivolvement of electrostatic channeling in guiding substrate transfer between two of the enzymes of 
the citric acid cycle. Similar effects may operate between other members of the citric acid metabolon. 



The association of enzymes involved in a particular 
metabolic pathway into a single multienzyme complex is a 
well-known phenomenon, with perhaps the most prominent 
example being the pyruvate dehydrogenase complex (Reed, 
1974) responsible for production of acetyl coenzyme A. The 
purported advantages of forming such multienzyme com- 
plexes seem intuitively reasonable (Welch, 1977; Srere, 
1985): the close proximity of enzymes responsible for 
catalyzing consecutive steps of a metabolic cycle may reduce 
the time required for substrate to diffuse between enzymes 
and prevent the substrate from escaping into solution where 
it might be sequestered by other enzynies for use in different 
metabolic pathways. This latter aspect may be particularly 
important for maintaining metabolic flux with low overall 
concentrations of metabolite (Welch, 1977). 

In principle, consecutive enzymes in a multienzyme 
complex might be close enough together that simple random 
diffusion would be sufficient to provide an efficient means 
of transferring substrate between active sites. Recent 
experimental and theoretical work suggests however that the 
efficiency of substrate transfer might be dramatically en- 
hanced by electrostatic effects. This idea stems largely from 
the solution of the crystal structure of the protozoan 
bifimctional enzyme dihydrofolate reductase- thy midy late 
synthase (DHFR-TS)' (Knighton et ai„ 1994). DHFR and 
TS catalyze sequential steps in die thymidylate cycle, with 
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the dihydrofolate (FH2) produced by the TS-catalyzed 
reaction subsequently acting as the substrate for the DHFR- 
catalyzed reaction. Experimentally the bifimctional form of 
the enzyme has been shown to have significant kinetic 
advantages over the related monofimctional enzymes in terms 
of a markedly reduced lag time for appearance of the final 
products of the coupled system (Meek et al., 1984). The 
crystal structure of DHFR-TS shows the entrances to the 
active sites to be at the enzyme surface and separated by 
around 40 A. Lying between the two active sites is a number 
of positively charged residues which were hypothesized to 
guide FH2, which has a net charge of -2e, fi-om the TS to 
the DHFR active site. Strong support for this idea has come 
from Brownian dynamics (BD) simulations of the process 
of substrate transfer (Elcock et ai, 1996). In simulations 
conducted at zero ionic strength the efficiency with which 
substrate was transferred from the TS to the DHFR active 
site was found to be almost 100%. Perhaps more impor- 
tantly, the simulations suggested that uncharged substrates, 
which experience no electrostatic guidance mechanism, are 
transferred with an efficiency of only around 5%: much too 
low to explain the experimental observation of efficient 
channeling. Further support for electrostatic channeling in 
DHFR-TS has been provided by recent experimental work 
showing that the efficiency of substrate transfer is dependent 
on ionic strength {Trujillo et aL, 1995, 1996); this result, 



' Abbreviations: CS. citrate synthase; MDH, malate dehydrogenase; 
PB, Poisson-Boltzmann; BD, Brownian dynamics; DHFR-TS. dihy- 
drofolate reductase-thymidylate synthase; NADH, nicotinamide ad- 
enine dinucleotide (reduced); AAT, aspartate aminotransferase. 
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which was also reflected in the simulations, is an expected 
characteristic of an electrostatic mechanism of transfer. 

Given the apparent importance of electrostatic effects 
suggested by both experimental and theoretical results for 
DHFR— TS, it is clearly of considerable interest to ask 
whether a similar mechanism operates elsewhere in nature. 
The focus of the present work is therefore on investigating 
the possible role of an electrostatic channeling mechanism 
in transferring substrate oxaloacetate between the citric acid 
cycle enzymes malate dehydrogenase (MDH) and citrate 
synthase (CS). Recent experimental work (Lindblah et al,, 
1994) has described the preparation and characterization of 
a fusion protein of the yeast forms of the two enzymes in 
which the C-termini of CS are connected by a short linker 
sequence (Gly-Ser-Gly) to the N-termini of MDH. Kinetics 
studies of the fusion protein resulted in two findings that 
suggest that oxaloacetate is efficiently channeled between 
MDH and CS. First, the lag (transient) time for formation 
of the final products (NADH and citrate) in the fusion protein 
was shorter than in the free enzymes; the lag time was also 
shorter than the theoretical estimate based on the ratio of 
Umax and Km for the CS part of the coupled system (Easterby, 
1973). Second, addition of aspartate aminotransferase 
(AAT), which competes with CS for the intermediate 
substrate oxaloacetate, had much less of an effect on the 
steady-state rate of the fusion protein than that of the free 
enzymes, an observation that indicates that oxaloacetate is 
largely sequestered from the bulk solution. 

From a simulation viewpoint the MDH— CS fusion protein 
is attractive since good structural information, a prerequisite 
for theoretical work, is available for each of the individual 
components of the system. While crystal structures of the 
yeast forms of the enzymes have not been solved, structures 
are available for the pig heart forms of both enzymes, and 
these are expected to be good models for the yeast enzymes. 
The problem of docking the two enzymes together in some 
reasonable fashion remains of course, but the use in the 
experimental work of a very short linker region (three amino 
acids) to connect the CS and MDH enzymes means that the 
orientational possibilities for such a fusion protein are 
drastically reduced. In fact, this sufficiently constrains the 
system that it is possible for a model of the fusion protein 
to be constructed by simply docking the CS and MDH 
structures so that the CS C-termini are adjacent to the MDH 
N-termini (Lindblah et al., 1994). 

A model of the fusion protein built in this way has the 
MDH and CS active sites separated by neariy 60 A, with 
the MDH active site facing away from CS. This is such a 
large distance that it would appear at first sight unlikely that 
substrate would be transferred between active sites with any 
great efficiency: certainly, it appears to have prompted Srere 
and co-workers to consider the possibility of a dimer of the 
fusion protein being responsible for efficient channeling 
(Lindblah et al., 1994). We report here, however, the results 
of BD simulations which suggest that, even in a monomer 
of the fusion protein, an electrostatically based channeling 
mechanism can be surprisingly efficient. 

MATERIALS AND METHODS 

All BD simulations were performed with the program 
UHBD (Madura et al., 1994, 1995). The crystal structures 
of MDH (Gleason et aL, 1994) and CS (Remington et al.. 
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Figure 1 : Ca trace of the MDH-CS fusion protein model. The 
CS dimer is shown in red, with the three C-terminal residues of 
each monomer colored purple. The MDH dimer is shown in blue, 
with the three N-terminal residues of each monomer colored green. 
This figure was prepared using the program QUANTA. 



1982) (both of which are dimers) were used as the starting 
point for generating a model of the fusion protein. The 
"open" form of CS was used as it has been suggested to 
represent a substrate entry/product release form of the 
enzyme (Remington et ai, 1982). Bound citrate molecules 
were removed from both the MDH . and CS structures. 
Hydrogens were added to each enzyme using the molecular 
simulation program CHARMM (Brooks et al., 1983). Ti- 
tratable residues were assumed to be in their usual proto- 
nation state at pH 7; i.e., the net charge of all Asp and Glu 
residues was set to — le and all Arg and Lys residues set to 
H-le; protonation states for histidines were assigned on the 
basis of their local environment, and all were in neutral 
forms. The total charge on the MDH dimer was therefore 
-f 6e; that on the CS dimer was — 2e. The two enzymes were 
docked so that the C-termini of the CS dimer were within 5 
A of the N-termini of the MDH dimer and so that structural 
overlap was kept to a minimum (Figure 1). A more refined 
method of docking involving energy minimization of the 
resulting structure could of course be envisaged, but the 
structural uncertainties involved are large enough at this stage 
that we feel that such an approach would not be justified; in 
any case, the results are not likely to be altered significantly. 

Prior to BD simulations being performed, the electrostatic 
potential around the fusion protein was obtained by finite- 
difference solution of the linear Poisson— Boltzmann equation 
(Honig & Nicholls, 1995) using a 160^ grid of spacing 1.2 
A. Boundary potentials were set assuming that each charge 
in the system is subject to Debye-Huckel screening. Atomic 
radii and charges were obtained from the CHARMM22 
parameter set (Mackerell et al., 1995). The internal dielectric 
of the protein was set to 4.0 while that of the solvent was 
set to 78.0, appropriate for water at 25 °C. The ionic strength 
was set to 0 mM in most simulations, although a range of 
values between 0 and 150 mM was used to investigate the 
ionic strength dependence of transfer efficiency. Oxalo- 
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Figure 2: MDH-CS system viewed fh)m a similar position to that in Figure 1 . The molecular surface of the protein is colored according 
to the electrostatic potential calculated at zero ionic strength. Blue areas on the molecular surface represent electrostatically positive regions; 
red, electrostatically negative regions (limits set at ifc 1 5 kT). The green area indicates the substrate starting point at the first MDH monomer's 
active site. The purple area marks the substrate reaction region at the first (most favored) CS active site. The second (less favored active 
site) CS active site is around the back of the figure. Substrate molecules reaching the purple surface are assumed to react at the CS active 
sites. The large sphere marks a radius 90 A firom the center of the protein: substrate molecules which reach any point on this surface are 
considered to have escaped to bulk solution. This figure was prepared using the program GRASP. 



acetate was modeled as a sphere of radius 2 A and charge 
— 2e; as in previous work (Elcock et aL, 1996) this does not 
represent a particularly realistic description of the substrate's 
atomic level structure, but electrostatic channeling effects, 
if present, are likely to be largely accounted for by the overall 
charge of the substrate. Simple spherical models of substrate 
have worked well in other simulations, for example in studies 
of the diffusion of acetylcholine to the active site of 
acetylcholinesterase (Antosiewicz et al., 1995). A series of 
BD simulations were performed starting from 1 180 different 
positions evenly distributed over a surface 12 A from the 
position of the CI atom of the bound citrate molecule in the 
active site of the first MDH monomer. These starting points 
are located adjacent to the flexible loop formed by residues 
77—87 known to play a role in substrate binding (Gleason 
et aL, 1994); these points therefore appear to be a reasonable 
model for the location of oxaloacetate emerging from the 
MDH active site. Substrate trajectories were simulated using 
the Ermak— McCammon algorithm (Ermak & McCammon, 
1978) with a time step of 0.05 ps. Trajectories were 
propagated until one of the following criteria were satis- 
fied: (1) the substrate "reacted" at the CS active sites; (2) 
the substrate reached a distance of 90 A firom the center of 
the fusion protein (at which point it was assumed to have 



escaped to bulk solution (see Figure 2)); (3) the trajectory 
exceeded 1 x 10^ steps. The last criterion was required 
because a few simulations remained trapped at the MDH 
active site where they would have stayed indefinitely; such 
trajectories were omitted from the analysis of results. A 
"reaction" was assumed to occur when the substrate came 
within 8 A of the position of the C3 atom of the bound citrate 
molecule observed at the active sites of the CS crystal 
structure (Remington et al,, 1982). An 8 A distance was 
selected since it gives a reaction region small enough that it 
provides a good definition of the active site area but 
sufficiently large that it is unaffected by changing the 
molecule's grid positioning; the results reported here are not 
sensitive to this choice of reaction distance. Substrate 
transfer efficiency was calculated as the percentage of 
trajectories which terminate successfully at one or other of 
the CS active sites. 

RESULTS 

In the simulations, the efficiency of substrate transfer from 
the MDH active site to either of the CS active sites is strongly 
dependent on the net charge of the substrate (Figure 3), 
dropping sharply as the charge is changed fi-om — 2e (that 
of oxaloacetate) to +2e. A more important result is that 
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Figure 3: Dependence of substrate transfer efficiency on substrate 
charge at zero ionic strength. Transfer efficiency is defined here 
as the percentage of substrate molecules leaving the MDH active 
site which reach a CS active site. 




Figure 4: Electrostatic isopotential contours around the MDH- 
CS fusion protein (calculated at zero ionic strength) viewed from 
the same position as in Figure 2. The blue surface connects all 
points having a potential of + 1.5 kT; the red surface connects all 
points having a potential of -1.5 kT. Note the presence of a 
continuous region of positive electrostatic potential covering the 
area where the MDH and CS dimers meet. 

the transfer efficiency for an uncharged substrate, which is 
not subject to any electrostatic channeling effects (and 
therefore behaves in exactly the same way as a charged 
substrate would with the enzyme's electrostatic field switched 
off), is only around 0.5 ± 0.4%, This is perhaps what one 
would expect, given the relative positioning of the active 
sites in the model of the fusion protein (Figure 2). In 
contrast, the transfer efficiency for a charge of — 2e is 43 it 
2%. This result, which does not seem at all obvious fi^om 
simple inspection of the model structure, represents an 
approximate 100- fold increase in efficiency over the simula- 
tions without electrostatic channeling effects. On the other 
hand, the result is less surprising when the electrostatic 
potential around the fusion protein is examined (Figure 4); 
isopotential contours marked at values of ±1.5 kT indicate 



the presence of a large area of positive (i.e., favorable) 
electrostatic potential extending over the region where MDH 
and CS meet. That the transfer efficiency does not exceed 
50% appears to be due to the fact that the region of favorable 
potential does not quite extend to the substrate starting 
positions, though the apparently unfavorable positioning of 
the MDH active site probably also plays a role. Once 
substrate molecules reach the region of favorable electrostatic 
potential however, they are likely to remain within it and, 
once trapped in this way, stand a good chance of ultimately 
diffusing to one or other of the CS active sites. 

Simulations performed to assess the relative contributions 
of the two enzymes in determining the efficiency of substrate 
transfer suggest that both play important roles, though CS 
dominates. When the electrostatic field due only to MDH 
is switched off (i.e., when all atomic charges on MDH are 
set to zero), the transfer efficiency is reduced by a factor of 
4 to 12 ± 2%; when the electrostatic field due only to CS is 
switched off, the transfer efficiency is reduced even more 
dramatically, to only 0.6 ± 0.5%. This apparent dominance 
of CS's electrostatic field was further emphasized by the 
results of simulations run in the opposite direction, i.e., 
starting at the CS active site and terminating at the MDH 
active sites. In the DHFR— TS system such "backward" 
simulations were of comparable efficiency to the "forward" 
results, suggesting that in DHFR— TS the role of electrostatics 
is more one of restricting the substrate's diffusion to the 
region between the active sites than one of actively guiding 
substrate from one active site to the other (Elcock et al.y 
1996). Results for the MDH-CS system are more difficult 
to interpret since all backward trajectories remained trapped 
indefinitely at the CS active site, making it impossible to 
calculate substrate transfer efficiencies. However, since this 
trapping results from strong electrostatic effects (no trapping 
was observed when the enzyme's electrostatic field was 
switched off), it is tempting, and perhaps not unreasonable, 
to interpret this result as implying that, in contrast to DHFR— 
TS, the electrostatic channeling in the MDH— CS system is 
indeed directional. 

Consistent with the above results suggesting the impor- 
tance of electrostatic effects, the transfer efficiency obtained 
from the simulations is found to be strongly dependent on 
the ionic strength of the solution (Figure 5); at 150 mM the 
efficiency drops to a value of 1 1 ± 2%. Again, this result 
is understandable in terms of the electrostatic potential 
(Figure 6); the increased shielding of the fusion protein's 
charged residues means that the regions of favorable 
electrostatic potential no longer extend out so far into solution 
so that the probability of substrate escape is increased. 

In the model that we have constructed substrate molecules 
start their trajectories at one end of the fiision protein (Figure 
2). Since it is possible for the substrate to diffuse from this 
position around either side of the protein, reactions are 
expected to be obtained at both of the CS active sites. 
However, as the substrate starting positions are slightly 
displaced toward one side of the protein, it is perhaps not 
surprising that an approximately 2-fold preference is obtained 
for reaction at one of the CS active sites (from here on, we 
use "first" active site to mean "more favored") (Figures 3 
and 5). More interesting is the fact that the distributions of 
reaction times (i.e., trajectory lengths) for reactions at the 
two CS active sites are very different (Figure 7). For the 
less-favored CS active site, the probability of a particular 
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Figure 5: Dependence of substrate transfer efficiency on ionic 
strength for a substrate charge of -2. 




Figure 6: Electrostatic isopotential contours around the MDH— 
CS fusion protein (calculated at 1 50 mM) viewed from the same 
position as in Figure 2. The blue surface connects all points having 
a potential of + 1.5 kT; the red surface connects all points having 
a potential of -1.5 kT. 

reaction time increases monotonically as a function of the 
reaction time itself In contrast, for the more favored CS 
active site, a distinct maximum in the distribution is observed 
at much shorter times, superimposed on a trend similar to 
that seen with the second CS active site (Figure 7). in other 
words, a considerable proportion of the reactive trajectories 
for the first CS active site reacts very rapidly: out of 438 
trajectories which terminate at the first CS active site, 113 
(26%) do so in less than 10 ns, whereas only 4 out of 227 
reactive trajectories (2%) at the second CS active site react 
within this time. The reasons for this dramatic difference 
are not immediately obvious, but it is interesting to note that 
no correlation appears between the substrate's starting 
position and either its eventual fate or its speed of reaction; 
any concern that the results might be very sensitive to the 
former aspect can therefore be ruled out. Instead, it appears 
likely that the bimodal distribution reflects real differences 
in substrate trajectories, though it is by no means clear of 
course whether these effects will ultimately be of much 
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Figure 7: Distribution of reaction times (i.e., trajectory lengths) 
for reaction at the two CS active sites obtained at zero ionic strength. 

relevance to overall channeling efficiency; further work will 
be required to address this problem more closely. 

DISCUSSION 

The simulation results reported here provide evidence that 
the experimental observation of efficient substrate channeling 
in the MDH— CS fusion protein is due to an electrostatic 
mechanism. Since electrostatic effects appear so important 
in the simulations, it is of course no surprise to find that if 
channeling is very efficient for a substrate charge of — 2e 
(as it is here), it is extremely inefficient for a substrate charge 
of H-2e. On the other hand, comparison with results obtained 
for a substrate charge of zero is more illuminating since it 
gives not only a measure of the increase in transfer efficiency 
that might result from electrostatic channeling but also an 
indication of whether electrostatic channeling is required at 
all in the first place. In the present case, the MDH and CS 
active sites are separated by around 60 A, and in the absence 
of electrostatic effects (i.e., when the substrate has a charge 
of zero) much less than 1% of substrate molecules leaving 
the MDH active site reach the CS active site. The simula- 
tions suggest therefore that efficient channeling of substrate 
does not result simply from the proximity of the CS and 
MDH active sites. When electrostatic effects are included 
however (i.e., when the substrate has a charge of — 2e), the 
transfer efficiency not only increases by around 2 orders of 
magnitude but also approaches a value that is more consistent 
with the experimental observation of efficient substrate 
channeling. 

Many of the results reported here for the MDH- CS system 
parallel those we obtained previously in our studies of 
channeling in DHFR-TS (Elcock et ai, 1996) and are 
obviously those to be expected of an electrostatic channeling 
mechanism. The ionic strength dependence of transfer 
efficiency obtained here appears more marked than in our 
study of the DHFR— TS system, but not surprisingly follows 
the same qualitative trend: increasing ionic strength tends 
to suppress any electrostatic channeling effect which might 
be present. The magnitude of this ionic strength effect is 
likely to vary from system to system and, in extreme cases 
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of very strong electrostatic effects, might even be almost 
completely insensitive to ionic strength within the physi- 
ological range. However, in the present case it seems that 
ionic strength changes should result in experimentally 
measurable changes in the kinetics of the MDH— CS protein. 

This aspect leads us to consider how contact between our 
theoretical fmdings and experimental results is to be made; 
it is important to stress, for example, that there is as yet no 
direct experimental evidence that transfer of oxaloacetate is 
electrostatically mediated. One problem with connecting the 
present work to experimental results is the use here of pig 
heart enzyme structures, when the experimental work 
involved the yeast forms of the enzymes. As stated in the 
introduction, this difference is unavoidable owing to the lack 
of structural information for the yeast enzymes. A related 
difficulty involves our overall structure of the fusion protein; 
in the absence of further structural information we have no 
way of knowing whether our structure represents a good or 
poor model of the real fusion protein. Before dismissing 
the results reported here as artifacts of our modeled structure 
however, it should be borne in mind that it places the MDH 
active site facing away fi-om CS and thus would not be 
expected to bias the results toward very efficient channeling. 
A more fundamental problem, and one of obvious impor- 
tance, is how to connect the transfer efficiencies obtained 
by simulation with lag times and other observables obtained 
by experimental kinetics studies. Preliminary results (un- 
published work) based on a reaction kinetics scheme suggest 
that the transfer efficiency of ~45% obtained in our 
simulations is probably sufficient to account for the chan- 
neling of oxaloacetate observed in the AAT competition 
experiments (Lindblah et ai., 1994). The development of a 
more complete method for connecting theoretical and 
experimental results, which will clearly be essential for 
example for interpreting ionic strength effects, will form the 
subject of further work. 

Irrespective of how directly our fmdings can be compared 
to experimental results, the 100-fold increase in transfer 
efficiency obtained by changing the substrate charge from 0 
to — 2e in our simulations argues strongly that electrostatics 
may be of importance for the MDH— CS system. It is 
interesting therefore to ask whether such a mechanism is also 
likely to operate elsewhere in the citric acid cycle, particularly 
in view of the fact that all of the cycle's substrates 
(intermediates) have a charge of — 2e or — 3e. The existence 
of a multienzyme complex (metabolon) formed by enzymes 
of the citric acid cycle has been known for some years 
(Barnes & Weitzman, 1986), though it remains poorly 
characterized (Robinson et al., 1987). From a simulation 
point of view, interesting extensions of this work would 
involve investigating the two enzymes bracketing MDH and 
CS in the citric acid cycle, namely, fumarase and aconitase. 
High-resolution structures are available for both enzymes, 
but further information would be required before further 
theoretical work would be justified: in the absence of 
additional structural constraints, the large number of possible 
relative orientations of the enzymes would preclude mean- 
ingful study. 

While the options for further theoretical study of chan- 
neling effects in multienzyme complexes are limited by the 
requirement for good structural information, experimental 
studies, such as kinetics measurements, suffer from no such 
limitation. However, in the absence of structural data, the 



interpretation of the experimental results may well be 
difficult. It is important to point out for example that in 
studies in which two enzymes, which would otherwise 
remain separate, are forced to associate [for example by the 
use of a covalent linker as in the work of Lindblah et al, 
(1994)] it may not be possible to decide with certainty 
whether the absence of a channeling effect is real or simply 
due to the fact that structural constraints force the enzymes 
to adopt a relative orientation unfavorable for channeling. 
By the same token, the efficient channeling obtained in the 
MDH— CS fusion protein might conceivably result fi-om the 
ens^ymes being forced to adopt a relative orientation which 
just happens to be favorable for channeling. In the citric 
acid cycle metabolon itself, the two enzymes will probably 
be oriented in a different way, which may either increase or 
decrease the efficiency of channeling. 

Finally it is worth pointing out that, even if accurate 
structural data are available, the simulation results reported 
here indicate that simple inspection of a structure may not 
be a good guide to estimating the efficiency of substrate 
channeling. If an electrostatic mechanism is important, as 
is suggested strongly by the present results, the spatial 
proximity of active sites need not correlate at all with the 
substrate transfer efficiency; what is likely to be more 
important is the presence of an unbroken electrostatically 
favorable region connecting the active sites. 
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A gene encoding a bifunctional homodimeric dihydrofolate 
reductase-thymidylate synthase (DHFR-TS) was con- 
structed by destroying the stop codon of Escherichia coii 
dihydrofolate reductase (DHFR) and joining the coding 
sequences of the monofunctional enzymes by a five amino 
acid linker. The protein was designed to mimic features of 
active site proximity and electrostatics in the protozoan 
DHFR-TSs which are believed to be important in channeling 
of the DHFR substrate, Hjfolate, to TS. The genetlcaUy engi- 
neered catalytically active homodimeric bifunctional 
DHFR-TS was expressed, purified and characterized. The 
component activities of the purified bifunctional enzyme had 
kinetic properties similar to those of the monofunctional TS 
and DHFR, but unlike the authentic bifunctional enzymes 
from protozoa this enzyme did not kinetically channel dihy- 
drofolate from DHFR to TS. 

Keywords: biftinctional DHFR-TS/engineered DHFR-TS/ 
linker/kinetic channeling/substrate channeling. 



Introduction 

Thymidylate synthase (TS; EC 2,1.1.45) and dihydrofolate 
reductase (DHFR; EC 1.5.1.3) catalyze sequential reactions in 
the de novo synthesis of thymidine-5' -monophosphate (dTMP) 
(Figure 1). TS catalyzes the reductive methylation of 2'- 
deoxyuridine-5 '-monophosphate (dUMP) by 5,10-methyiene- 
5,6,7,8-tetrahydrofolate (CH2H4folate) to produce dTMP and 
7,8-dihydrofolate (Hifolate), and DHFR catalyzes the sub- 
sequent NADPH-dependent reduction of H2folate to generate 
5,6,7,8-tetrahydrofolate (H4folate). Both TS and DHFR exist 
as distinct, monofunctional enzymes in sources as varied 
as bacteriophage, bacteria, fungi, viruses and mammals 
(Ackermann and Potter, 1949; Blakley, 1984). In most genera, 
TS is a dimer of identical subunits of about 35 kDa (Carreras 
and Santi, 1995) whereas monofunctional DHFR is a monomer 
of about 20 kDa (Blakley, 1984). In contrast, TS and DHFR 
exist on the same polypeptide chain in protozoa (Ivanetich 
and Santi, 1990) and higher plants (Lazar ei qL, 1993). The 
bifunctional protein is a dimer of identical subunits with an 
N-terminal DHFR domain connected to a TS domain by a 
linker of variable length depending on the source of the enzyme. 
Substrate channeling has been proposed as a biological 
advantage of bifunctional DHFR-TSs over their monofunc- 
tional counterparts (Meek et ai, 1985). Kinetic evidence for 
substrate channeling of the H2folate produced by TS to the 
DHFR active site has been reported for the Leishmania major 
and Toxoplasma gondii DHFR-TSs (Meek et at., 1985; Trujillo 
et aLy 1996), but the mechanism of channeling remains 
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unknown. It is possible that the channeling of H2folate in the 
bifunctional proteins occurs because of the proximity of the 
TS and DHFR active sites that results from linkage of the 
domains. In this case, DHFR would simply be exposed to the 
higher local concentration of TS-generated H2folate before it 
dilutes into bulk solvent. However, from the recently reported 
crystal structure of Lmajor TS-DHFR, it has been suggested 
that the TS and DHFR sites may not be sufficiently close to 
permit channeling by proximity (Knighton et al, 1994). 
Another proposed mechanism for chaxmeling is based on the 
surface charge distribution of Lmajor DHFR-TS (Knighton 
e( ai, 1994). The protein has surfaces of positive electrostatic 
potential that lie between and within the TS and DHFR domain 
binding sites for H2folate, and the negatively charged H2folate 
could be shuttled directly between the two active sites by way 
of electrostatic interactions. Tlie crystal structures of the 
Escherichia colt DHFR and TS monofunctional enzymes reveal 
similar regions of positive charge around the active sites of 
the monofunctional enzymes (Stroud, 1994), Linking the 
monofunctional Exoll DHFR and TS to produce a bifunctional 
protein provides an approach to testing the aforementioned 
hypotheses for channeling. Such a protein could provide high 
H2folate concentrations near the DHFR site and also simulate 
the electrostatic environment surroimding the H2 folate binding 
sites, although it would lack some of the electrostatic properties 
of protozoan bifunctional DHFR-TSs which may contribute 
to the proposed channeling mechanism (Knighton et aL, 1994). 
Furthermore, it has been suggested that at the high local 
enzyme concentrations present in vivo (Srere, 1987), the 
monofunctional enzymes might associate (Kawai and Hillcoat, 
1974; Mathews et ai, 1988; Moen et al, 1988) and also be 
capable of channeling H2folate (Stroud, 1994). This hypothesis 
could also be tested by linkage of the two monofunctional 
enzymes into a bifunctional protein. 

Recently the coding sequences of the E.coli DHFR and TS 
were linked in an expression vector in an attempt to produce 
the bifunctional DHFR-TS homodimer (Iwakura and Kokubu, 
1995). Upon expression, TS homodimer and TS^DHFR-TS) 
heterodimer were isolated without the desired DHFR-TS 
homodimer. In the present paper, we describe the construction 
of a gene for a bifunctional Exoli DHFR-TS and its expression 
and purification as a catalytically active homodimer similar to 
the natural protozoan enzymes. We present the kinetics of 
H2folate transfer between the TS and DHFR domains of 
the bifunctional enzyme and report that there is no kinetic 
channeling of H2folate, 

Materials and methods 

Materials 

Escherichia coii strain DH5a, obtained from BRL, was used 
for general manipulation of recombinant plasmids. Escherichia 
coii strain Thy" x2913RecA (AthyA 572, recASS) was used 
for expression of the engineered bifunctional DHFR-TS (Cli- 
mie et al, 1992). pThyA, which contains the Exoli TS gene 
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Fig. 1. The reactions catalyzed by thymidylate synthase (TS) and 
dihydro folate reductase (DHFR). 

in Biuescript KS (Maley aiid Maley, 1988), was a gift from 
F.Maley. The expression vector pfrc 99 A, which contains the 
lac I** gene and the trc promoter, was purchased from Pharma- 
cia, General methods for ligation and transformation of E.coli 
were performed as previously described (Sambrook et al, 
1989). Restriction endonucleases, T4 ligase and vent DNA 
polymerase were from New England Biolabs, Life Technolo- 
gies, Inc. or Boehringer Mannheim. Oligonucleotides were 
synthesized at the Biomolecular Resource Center at the Univer- 
sity of California, San Francisco, using cyanoethyl phos- 
phoramidite chemistry. Automated DNA sequencing was 
performed at the Biomolecular Resource Center, University of 
California, San Francisco, using the Taq dye terminator cycle 
sequencing chemistry. CH2H4folate (Bruice and Santi, 1982), 
enantiomerically pure H2 folate (Roos, 1993) and methotrexate 
(MTX)-Sepharose CL-6B (~l ^mol MTX/ml resin) (Meek 
et al, 1985) were prepared as described. Purified E.coli DHFR 
and Exoli TS were gifts of C.R.Matthews and F.Maley, 
respectively. Protein concentrations were determined by the 
Bradford method using bovine serum albumin (BSA) as a 
standard (Bradford, 1976). SDS-PAGE was performed as 
previously described (Read and Northcote, 1981). 

Structure alignments 

The UCSF MidasPlus software system from the Computer 
Graphics Laboratory, University of California, San Francisco, 
was used to obtain alignments of the crystal structures of 
E.coli TS with 5-fluoro-2'-deoxyuridine 5 '-monophosphate 
(FdUMP) and 1 0-propargyl-5,8-dideazafolate (CBS 7 1 7) bound 
(Fauman et ai, 1994) and E.coli DHFR with folate and 
NADP-^ bound (Bystroff et ai, 1990) to the structure of 
L.major DHFR-TS with FdUMP and CB3717 bound to the 
TS domains and MTX and NADPH to the DHFR domains 
(Knighton et al., 1994). The structures were matched by means 
of a least-squares fit of the coordinates of the backbone a- 
carbons of three conserved active site residues and one 
additional conserved residue. The backbone a-carbons of the 
following residues were used for the TS alignment: C146, 
1264, E58 and L38 for Exoli and C400, 1520, E292 and L271 
for L.major, The backbone a-carbons of the following residues 
were used for the DHFR alignment: G15, W22, T35 and F125 
for Exoli and G40, W47, T60 and F193 for L.major. The 
distance between the DHFR C-terminus and the TS N-terminus 
of the matched monofunctional enzymes was then determined. 

Construction of expression vector pTDDHFR-TS 
The 5' end of the E.coli TS gene was modified by PCR using 
pThyA as a template. Two synthetic primers PI (5'- 
ACCG CTCGAG CGTCGTGGTGCGAGCTCTGGT GAA AA 
ACAGTATTTAGAACTG-3') and P2 (5'-CTGGTTCAG- 
TACCGTAGT-3') were used for the PCR amplification step. 
The PI primer introduced several features: the C-terminus of 
the coding sequence of DHFR with the stop codon removed, 
a linker region and the N-terminus of the coding sequence of 
TS with the first methionine changed to a glutamate (underlined 
and bold). In addition, an Xho\ site (underlined) for cloning 



and a Sad site (bold) within the linker were introduced. P2 
primed in the reverse direction on the TS coding sequence 
downstream from the 55jHII site. The resulting DNA fi^gment 
was digested with Xliol and fi^^HII and cloned into pThyA to 
yield pTDI, which was checked by restriction analysis and 
DNA sequencing of the amplified region. 

The Exoli DHFR coding sequence was amplified by PCR 
using DH5a chromosomal DNA as a template. Two synthetic 
primers P3 (5'-CCCCGGTACCTCATGATCAGTCTGATT- 
GCG-3') and P4 (S'-CCGGCTCGAGAATCTCAAAGCAAT- 
AGCT-3') were used. The P3 primer introduced a Kpn\ site 
(underiined) and a BspWl site (bold) at the beginning of the 
DHFR coding sequence. P4, which includes a unique Xhol 
site (underlined) present in the DHFR coding sequence, primed 
in the reverse direction two amino acids before the C-terminus 
of the DHFR coding sequence. The resulting DNA fragment 
was digested with Kpnl and Xhol and then cloned into pTDl 
to give pTD2, which was checked by restriction analysis and 
DNA sequencing of the amplified region. pTDI and pTD2, 
both of which are in Biuescript KS, do not contain ribosome 
binding sites. In order to achieve controlled expression of the 
biftinctional gene, it was necessary to transfer the construction 
into another vector such as ^Trc 99A. Digestion of pTD2 with 
BspHl and BomWl followed by ligation of the resulting 
fragment with p7Vc 99 A that had been restricted with yVcoI 
and BamV{\ yielded the final construction, pTDDHFR-TS. 

Expression and purification of Exoli DHFR-TS 
pTD DHFR-TS was transformed into E.coli strain Thy" x2913. 
A suspension of cells from 3 1 IPTG-induced culture was lysed 
using a French press. High molecular weight nucleic acids 
were precipitated by the addition of streptomycin sulfate. The 
resulting supernatant was loaded at 4°C on a Q-Sepharose 
column (30 ml bed volume), which had been pre-equilibrated 
with buffer A (20 mM Tris-HCl pH 7.4, 1 niM EDTA and 
1 mM DTT). The column was washed extensively, with buffer 
A containing 0,3 M KCl. The protein was eluted with a 300 
ml linear gradient of 0.3 to 1 M KCl in buffer A. Escherichia 
coli DHFR-TS eluted at approximately 0.7 M KCl. KCl was 
added to the pooled fractions to a final concentration of 
1,2 M KCl. The pool was loaded at 4°C on a phenyl- 
Sepharose column (30 ml bed volume) that had been previously 
equilibrated with buffer A containing 1.2 M KCl, The protein 
did not bind well to the column and it was recovered in an 
initial wash of the column with buffer A containing 1.2 M 
KCl. The protein was concentrated using Centriprep-30 from 
Amicon to a final volume of 10 ml. 

Enzyme assays 

Thymidylate synthase activity was assayed spectrophoto- 
metrically by monitoring formation of H2folate at 340 nm 
(£340 = 6400 M-' cm-') in 50 mM TES (pH 7.8), 0.1 mM 
CH2H4folate, 0.1 mM dUMP, 75 mM 2-mercaptoethanoI, 
1 mM EDTA and 5 mM HjCO (Wahba and Friedkin, 1961), 
DHFR activity was monitored by measuring the decrease in 
absorbance at 340 nm resulting from consumption of H2folate 
and NADPH (£340 = 12 300 M"' cm"') in 50 mM TES 
(pH7.0), 0.1 mM Hjfolate, 0.1 mM NADPH, 75 mM 2- 
mercaptoethanol, I mM EDTA and 1 mg/ml BSA (HiUcoat 
et al., 1967). For both TS and DHFR, one unit of enzyme 
activity is defined as the amount of enzyme that produces 
1 ^mol of product per minute at 25X. 

For determination of the steady-state kinetic parameters, the 
concentration of one substrate was varied between approxi- 
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mately 0.5 and 10 K^^, with the other substrate fixed at 100 ^M. 
The Kj„ value for H2 folate was determined using a 10 cm 
cuvette. Kinetic constants were calculated from a non-linear 
square fit to the Michael is~Menten equation using Kaleida- 
graph (Abelbeck Software) on a Macintosh computer. 

Inhibition of DHFR by pyrimethamine (Pyr) was measured 
by varying the concentration of inhibitor in a reaction mixture 
containing 4.3 nM enzyme and 10 ^iM H2folate. The Ki was 
calculated using the equation for competitive inhibition from 
Segel (Segel, 1975). 

Coupled assay 

The DHFR-TS coupled system was assayed by monitoring 
the decrease in NADPH absorbance at 340 nm in the absence 
of added H2folate (Meek et ai, 1985). The assay mixture 
containing 28 ^iM CHiajfolate and 30 )xM NADPH in 50 mM 
TES (pH 7.8), 75 mM 2-mercaptoethanol, 1 mM EDTA, 5 mM 
H2CO, 1 mg/ml BSA and variable amounts of enzyme was 
incubated at room temperature until no change in absorbance 
at 340 nm was detectable. In this assay the formaldehyde 
reacts with H4folate to regenerate CH2H4folate, which makes 
the assay cyclic. The reaction was initiated by addition of 25 
|il of 8 mM dUMP (final concentration 200 |aM), and the 
decrease in absorbance was monitored at 340 nm (8340 = 5900 
M~' cm~'). 

In addition, individual TS and DHFR rates of the bifunctional 
enzyme were determined under conditions of the coupled 
assay, with the following modifications: for DHFR, dUMP 
was omitted and 100 |iM H^folate was included; for TS, 
NADPH was omitted. 

Analytical gel filtration 

Gel filtration of E.coli DHFR-TS was performed on a Superose 
12HR 10/30 column. One hundred |il of the purified enzyme 
(1 mg/ml) in buffer A containing 1.2 M KCl was injected, 
and the column was eluted with buffer A at a flow rate of 0.5 
ml/min. The column was calibrated with molecular weight 
standards of 43 to 2000 kDa. 

Results 

Structure alignments 

The mono functional crystal structures of Exoli TS (Fauman 
et ai, 1994) and Exoli DHFR (Bolin et ai, 1982) were 
matched to that of Lmajor DHFR-TS (Knighton et al., 1994) 
using the UCSF MidasPlus program. The alignment of the TS 
domains yielded a root mean square error of 0.42 A and that 
of the DHFR domains gave an error of 0.67 A. A similar 
alignment of the monofunctionai E.coli enzymes with the 
Lmajor protein has been performed by Stroud (see Figure 2 
of Stroud, 1994) who compared the electrostatic potential 
obtained for the E.coli enzymes to that obtained by Knighton for 
the Lmajor DHFR-TS (see Figure 6 of Knighton et ai, 1994). 

The interface of the matched E.coli structures is shown in 
Figure 2 along with that of the Lmajor biftinctional enzyme. 
The fit of the two surfaces was achieved by choosing different 
side chain rotomers of some interface residues. In this model, 
there are no steric clashes or polarity problems apparent 
between sidechains of the two monofunctionai enzymes. Two 
sets of surface sidechains of opposite polarity are close enough 
together that it would be possible for them adopt conformations 
that would permit charge interactions across the interface 
(DHFR D157 with TS R225 and DHFR KI09 with TS D253). 

The modelled distance between the monoftmctional DHFR 



C-terminus and the monofunctionai TS N-terminus was deter- 
mined to be about 10 A. A linker composed of five amino 
acids (GASSG) was used to connect the monoftmctional E.coli 
DHFR and TS enzymes (Figure 3). The distance between 
the C- and N-terminus of the linker in its fully extended 
conformation was about 18 A, which is sufficiently long to 
allow the engineered E.coli DHFR-TS bifunctional enzyme to 
adopt the structure of the Lmajor bifimctional enzyme. 

Construction of expression vector pTDDHFR-TS 
The E.coli bifunctional DHFR-TS gene was successftilly 
cloned under the control of the inducible trc promoter in p7>r 
99 A. The gene includes the coding sequence for E.coli DHFR 
with the stop codon removed, the coding sequence for the 
linker region and the coding sequence of E.coli TS with the 
TS starting methionine changed to a glutamate. This amino 
acid was selected because all TSs have either a starting 
methionine or an acidic residue at this position. 

Expression and purification of E.coli DHFR-TS 
Crude extracts of IPTG induced cultures containing 
pTDDHFR-TS were analyzed by SDS-PAGE. A protein band 
representing -10% of the total soluble protein was detected 
which migrated with a molecular mass of 48 kDa, the predicted 
size of the bifuncfional protein. The soluble crude extract 
showed TS and DHFR activities of 0.08 U/mg and 0, 1 U/mg, 
respectively, whereas no TS activity was detectable in the 
control crude extract from cells harboring pTrc 9 9 A. 

Initial attempts to purify the bifunctional protein using a 
MTX-Sepharose column were unsuccessful. The protein bound 
very tightly and could only be eluted by 0.01 M NaOH, which 
destroyed its activity. Instead, the Exoli DHFR-TS was 
purified by treatment of the crude extract with streptomycin 
sulfate followed by sequential chromatography on Q- 
Sepharose and phenyl-Sepharose columns. A total of 16 
milligrams of apparently homogeneous protein was obtained 
(Figure 4), with TS and DHFR activities of l.l U/mg and 2.5 
U/mg, respectively. 

Characterization of Exoli DHFR-TS 

The purified enzyme migrated as a single 48 kDa band on 
SDS-PAGE (Figure 4). It eluted from Superose I2HR 10/30 
with an apparent molecular mass of -100 kDa, indicating that 
the native protein was a dimer of identical subunits. 

Kinetic parameters for DHFR and TS were determined for 
the biftanctional protein (Table I). The value for DHFR 
was 2.0 ^iM for H2 folate and the ^cat was 32 s"'. The K,^ value 
for TS was 4.3 \iM for CH2H4folate and the fr^at was 1.0 s~'. 
The IC50 value for Pyr was determined in a reaction containing 
10 \xM H2folate. The K, was calculated to be 47 nM. These 
parameters are similar to those previously reported for the 
E.coli monofunctionai enzymes (Stone and Morrison, 1986; 
Flowell et a/., 1990; Zapf al., 1993). 

Coupled assay 

A previously described coupled assay system (Meek et ai, 
1985) was used to monitor the DHFR-catalyzed reduction of 
the H? folate formed by TS. The DHFR-catalyzed reaction 
converts NADPH to NADP**", and the resulting decrease in 
the concentration of NADPH can be monitored at 340 nm. 
Under steady state conditions, Equation I describes the time 
course of NADP"^ formation (Easterby, 1973; Rudolph et ai, 
1979; Meek et al, 1985). 

[NADP^] - V,/ + (v,/v2)/r,(e-v/^- - I) (1) 
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Fig. 2. (a) Stereoscopic image of the interface of E.coli DHFR and TS. DHFR is shown in green, and TS is shown in purple. The C-tcmiinus of DHFR and 
the N-terminus of TS are labeled. The structures of the monofunctional enzymes were aligned as described in the text. The structures are as described (Bolin 
et aL, 1982; Fauman et ai, 1994) except that different sidechain rotomers of some surface residues have been chosen to avoid steric clash (DHFR Q108, 
D87, El 34 and R158 and TS E250 and E223). (b) Stereoscopic image of the interface of l.major DHFR and TS shown for comparison. 



DHFR domain 




TS domain 



Fig. 3. Engineered bifunctional E.coli DHFR-TS witli linker (empty boxes). 
The first methionine of wild-type TS was changed to a glutamate (shaded 
box). 

In this equation, is the for Hz*""'^**^' v, is the rate of TS 
under conditions of the coupled assay, and is the rate of 
DHFR under near-saturating concentrations of substrates. In a 
plot of [NADP'*'] versus time, extrapolation of the linear 
portion of the curve to the horizontal axis provides the lag 
time (KJvj) that precedes attainment of the steady-state 
concentration of H2folate and extrapolation to the vertical axis 
provides the steady-state concentration {v\KJv2) of H2folate. 
In a coupled assay in which the two enzymes do not channel 




1 2 3 

Fig. 4. SDS-PAGE gci showing purification of the recombinant Exoli 
DHFR-TS. Lane I, crude extract; lane 2. Q-Scpharosc pool; lane 3, 
phenyl-Scpharose pool. 
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Table I. Summary of DHFR and TS kinetic parameters 




DHFR activity 






TS activity 




Enzyme 




(HM) for 


Ki (nM) 




(jiM) for 






Hjfolate 


for Pyr 




CH2H4folate 


DHFR-TS 


32 


2.0 


47 


1.0 


4.3 


Wild-type DHFR 


29° 




10'' 


NA<= 


NA'^ 


Wild-type TS 








3.7*1 


19d 



^Howell etaL, 1990. 
*'Stone and Morrison, 1986. 
*^NA, not applicable, 
^Zapf ef o/., 1993. 



Table II. Coupled assay using a mixture of monofunctional E,coli TS and DHFR and the bifunctional E.coli DHFR-TS homodimer 



Rate (|jM/min) Lag time (s) Steady-state [H2 folate] (pM) 



TS DHFR Calc. Exp. Calc. Exp. 



Bifunctional O.Il 0.55 218 217 0.40 0.43 

DHFR-TS 

Monofunctional 0.53 1.43 126 150 1.1 1 1.14 

DHFR + TS 



the intermediate substrate, the experimentally determined lag 
time and steady-state H2 folate concentration will match those 
values calculated using Equation 1, while in an assay in which 
the enzymes do channel, both values will be reduced from 
those calculated using Equation I. 

The coupled assay was first performed using a mixture of 
the monofunctional E.coli enzymes, such that the TS rate was 
0.53 ^iM/min and the DHFR rate was 1.43 ^iM/min (Table II). 
The experimentally determined values of 150 s for the lag 
time and 1 . 14 ^iM for the steady-state concentration of H2folate 
are similar to those calculated using Equation 1 (126 s and 
1. 11 \iM, respectively), confirming the validity of the equation 
under our experimental conditions for a non-channeling system. 
As a control for a channeling system, the bifunctional DHFR- 
TS from Tgoncfii (Trujillo et o/., 1996) was used in the same 
assay; in that case the lag time was below our detection limit 
of about 20 s (data not shown) indicating that H2folate channels 



between the enzyme's TS and DHFR active sites. The coupled 
assay was then carried out with the E.coli bifunctional DHFR- 
TS, using a preparation of the enzyme in which the TS rate 
was 0.1 1 ^M/min and the DHFR rate was 0.55 ^M/min (Figure 
5). In this case, both the experimentally determined and the 
calculated lag times were about 220 s, suggesting that there 
was no kinetic channeling of H^folate from the active site of 
TS to that of DHFR. The experimentally observed steady-state 
concentration of H2folate of 0,43 |iM was essentially the same 
as the calculated value of 0.40 |xM, confirming that the ftision 
of the DHFR and TS enzymes does not affect the rate of 
Hifolate transfer between the two active sites. 

Discussion 

It has been shown that both the L. major and Tgondii biftinc- 
tional DHFR-TSs channel the intermediate, H2folate, between 
the TS and DHFR active sites of the enzymes (Meek et ai. 
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1985; Trujillo et ai, 1996). In the crystal structure of the 
Lmajor protein, the active sites within a subunit are separated 
by a distance of about 40 A and those on opposite subunits 
are separated by about 70 A (Knighton et ai, 1994). Knighton 
and colleagues suggested that both distances appear to be too 
long to allow channeling of the intermediate on the basis of 
proximity. They proposed that a possible driving force for 
charmeling of H2folate could be electrostatic interactions 
between the negatively charged intermediate and the numerous 
positively charged patches on the surface of the protein. Stroud 
then noted that alignment of the structures of mono functional 
E.coli TS and DHFR with that of the Lmajor biftinctional 
enzyme revealed similar areas of positive potential (Stroud, 
1994). He suggested that E.coli TS and DHFR might be 
capable of channeling H^folate at the high local enzyme 
concentrations present in vivo either through the proposed 
electrostatic patliway or by way of a hinging mechanism which 
would allow the active sites to approach each other. An 
approach to testing this hypothesis is to join E.coli DHFR and 
TS via an appropriate linker and to determine whether the 
resulting engineered bifiinctional protein channels H2folate. 

In a recent study by Iwakura and Kokubu, the E.coli 
monofunctional DHFR and TS enzymes were joined by a short 
peptide linker which included a ribosome binding site (Iwakura 
and Kokubu, 1995). Because the TS gene was preceded by a 
ribosome binding site and because it contained a start codon, 
they obtained a mixture of products that included heterodimeric 
biftinctional TS-(DHFR-TS), homodimeric monofunctional 
TS and an intractable aggregated form of biftinctional DHFR- 
TS. Studies of the heterodimeric TS-(DHFR-TS) did not 
permit definitive conclusions to be drawn regarding kinetic 
channeling between the TS and DHFR sites. 

In the present study, a fiiUy active homodimeric biftinctional 
E.coli DHFR-TS was constructed by connecting the monoftmc- 
tional E.coli enzymes via a five amino acid long linker. The 
linker was designed to allow the engineered E.coli bifunctional 
DHFR-TS the possibility of adopting the conformation of the 
L.major biftinctional enzyme. The 18 A length of the ftilly 
extended linker was more than sufficient to span the necessary 
modeled 10 A distance. In addition, the linker was made 
flexible by the presence of glycines and soluble by the 
presence of serines (Robinson and Sauer, 1996). Expression 
and purification provided the desired homodimeric biftinctional 
protein with a subunit molecular weight of about 48 kDa and 
both TS and DHFR activities. Steady-state kinetic parameters 
of the individual enzymes of the bifunctional protein were 
similar to those observed for the wild-type monomeric E.coli 
DHFR and the dimeric E.coli TS. 

Kinetic channeling in the E.coli biftinctional DHFR-TS 
coupled system was studied by examination of the time 
required for the appearance of NADP'*' and calculation of the 
steady state concentration of H2 folate. For a coupled system 
in which there is no channeling, the lag time and the steady state 
concentration of intermediate can be predicted by Equation 1, 
while in a system that channels the experimentally determined 
values of both parameters will be reduced from those predicted. 
This approach has been used previously to demonstrate kinetic 
channeling of H2 folate for both Lmajor and T.gondii DHFR- 
TSs (Meek et ai, 1985; Trujillo et ai, 1996). The experiment- 
ally determined lag time and the steady-state concentration of 
H2 folate for the engineered E.coli TS-DHFR were similar to 
the values predicted for a mixture of monofiinctional DHFR 
and TS with the same rates, indicating that the biftinctional 



homodimeric DHFR-TS does not channel H2folate from the 
active site of TS to that of DHFR. 

Because our engineered bifiinctional E.coli DHFR-TS does 
not channel H2folate, we conclude that the local concentration 
of H2folate obtained by linking the monofunctional E.coli 
enzymes is not sufficient to promote substrate channeling. 
Although the E.coli DHFR-TS may possess some of the 
features of the electrostatic environment of the Lmajor enzyme, 
our results regarding the role of electrostatic facilitation of 
H2folate transfer between the active sites are inconclusive 
since some features which may be important in electrostatic- 
mediated channeling are absent in our protein. For example, 
a highly basic loop present in the DHFR domain of the Lmajor 
DHFR-TS and absent in our E.coli DHFR-TS has been 
proposed to play a role in the channeling of H2folate (Elcock 
et al., 1996). Finally, our finding that the conversion of the 
E.coli monoftinctional DHFR and TS to a biftinctional DHFR- 
TS does not result in kinetic channeling contradicts proposals 
that high concentrations of the monofunctional enzymes might 
facilitate substrate channeling in vivo. 
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Kinetics of the coupled reaction catalysed by a fusion protein of yeast 
mitochondrial malate dehydrogenase and citrate synthase 

Henrik Pettersson\ Peter Olsson^ Leif Bulow^ and Gosta Pettersson^ 

^Avdelningen fur Biokemi and ^Avdelningen fur Tilldmpad Biokemi, Kemicentrum, Lunds Vniversitet, Lund, Sweden 

The mechanistic implications of the kinetic behaviour of a fusion protein of mitochondrial malate dehydrogenase 
and citrate synthase have been reanalysed in view of predictions based on experimentally determined kinetic 
parameter values for the dehydrogenase and synthase activities of the protein. The results show that the time- 
course of citrate formation from malate in the coupled reaction catalysed by the fusion protein can be most 
satisfactorily accounted for in terms of a free-diffusion mechanism when consideration is taken to the inhibitory 
effects of NADH and oxaloacetate on the malate dehydrogenase activity. The effect of aspartate aminotransferase 
on the coupled reaction is likewise fully consistent with that expected for a free-diffusion mechanism. It is 
concluded that no tenable kinetic evidence is available to support the proposal that the fusion protein catalyses 
citrate formation from malate by a mechanism involving channelling of the intermediate oxaloacetate. 

Keywords: fusion protein; channelling; malate dehydrogenase; citrate synthase. 



The biological significance of a channelled transfer of 
metabolites between sequentially operating enzymes in meta- 
bolic pathways has been a matter of great dispute [1-3]. 
Several authors have argued that metabolite channelling is 
likely to occur due to 'proximity effects' once two sequential 
enzymes are brought close together, e.g. through the dynamic 
or static formation of a bienzyme complex. Such complex 
formation can be anticipated to decrease the average distance 
between the metabolite-producing and metabolite-consuming 
enzymic sites, and this has been envisaged to offer catalytic 
advantages in the form of enhanced steady-state and/or 
transient reaction rates [4-6]. 

The significance of proximity effects in enzymic metabolite 
transfer processes has been experimentally probed by kinetic 
studies of sequential enzymes that have been covalently fused 
to each others by chemical modification or genetic engineering 
[7-10]. Particular attention has been paid to fusion proteins 
containing the sequential enzymes malate dehydrogenase and 
citrate synthase in the metabolically central citric acid cycle. 
A fusion protein of the yeast form of these enzymes was 
biosynthesized and examined by Lindbladh et aL [11], who 
found that the transient lag time for product formation in the 
coupled reaction catalysed by the two enzymes was shorter for 
the fusion protein than for the free enzymes. This was taken to 
indicate that there is a channelled transfer of the intermediate 
oxaloacetate in the fusion protein, an interpretation that 
received support by the observation that inhibition of the 
coupled reaction by aspartate aminotransferase (which com- 
petes with citrate synthase for oxaloacetate) was less with the 
fusion protein than with the free enzymes. Similar results have 
later been obtained with a fusion protein of the porcine form of 
the enzymes [12]. 
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Enzymes: malate dehydrogenase (EC 1.1.1.37); citrate synthase 
(EC 4.1.3.7); aspartate aminotransferase (EC 2.6.1.1). 
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Unfortunately, kinetic parameter values for the malate 
dehydrogenase and citrate synthase activities of the fusion 
proteins differ somewhat from those of the free enzymes 
[11,12], This leaves some uncertainty as to whether the 
different kinetic behaviour of the fusion protein and the free 
enzymes may reflect the observed differences in kinetic 
parameter values rather than substrate channelling. To eliminate 
this uncertainty, we have now performed a more detailed 
kinetic study of the coupled reaction catalysed by the fusion 
protein containing the yeast form of the enzymes. The results 
do not support the previous interpretation of the reaction 
kinetics in terms of substrate channelling, but provide clear 
evidence that the kinetic behaviour of the fusion protein is fully 
consistent with a free-diffusion mechanism of oxaloacetate 
transfer between the malate dehydrogenase and citrate synthase 
moieties of the protein. 

EXPERIMENTAL PROCEDURES 

Materials 

Porcine heart mitochondrial malate dehydrogenase (MDH) and 
aspartate aminotransferase (AAT) were purchased from Sigma 
Chemicals Co. Crystals of the enzyme suspensions were 
collected by centrifugation and dissolved in 40 mM potassium 
phosphate buffer, pH 8.1, containing 0.1 mM EDTA and 
0.35 mM 2-mercaptoethanesulfonic acid. Immediately before 
use, enzyme solutions were passed through a HiTrap desalting 
column, pre-equilibrated with the above phosphate buffer. 
Active-site concentrations of the enzymes were determined 
spectrophotometrically at 280 nm using an absorption coeffi- 
cient of 19.8 mM~*.cm~* for MDH [13] and 65.8 mM"* cm'* 
for AAT [14]. 

The fusion protein of yeast citrate synthase and yeast 
mitochondrial malate dehydrogenase was prepared according to 
the protocol elaborated by Srere and coworkers, as detailed for 
their construction and purification of the fusion protein of the 
porcine enzymes [12]. Escherichia coli strain BL21, harbouring 
plasmid pODC29-Citl-3-Mdhl containing the gene for the 
fusion protein, was kindly provided by them. The purified 
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fusion protein obtained in the final gel filtration step was 
precipitated with 70% ammonium sulfate and stored at -20 ''C. 
Protein concentrations during purification were determined by 
the Bradford method [15]. The concentration of the fusion 
protein was estimated from its steady-state MDH activity in 
the direction of malate formation, as determined in 50 mM 
Tris/HCl buffer (pH 7.5) using 0.1 mM oxaloacetate and 
0.1 mM NADH. 1 unit of the fusion protein is defined as 
the amount of protein liberating 1 |xmol NAD"*" per min under 
such conditions. 

NAD"*" (grade I) was obtained from Boehringer Mannheim. 
Malate, oxaloacetate, acetyl-CoA, glutamate, aspartate, and 
dithionitrobenzoate were from Sigma Chemical Co. 

Methods 

All kinetic experiments were carried out at 25 °C in 40 mM 
potassium phosphate buffer, pH 8.1. The equilibrium constant 
for the oxidation of malate by NAD"*" was determined by 
reacting 10 mM malate with 4 mM NAD"*" in the presence of 
about 0.6 jjiM malate dehydrogenase. The amount of NADH 
formed after equilibration of the reaction solution was 
calculated from the 340 nm absorbance changes observed, 
using an absorption coefficient of 6200 M~* cm~*. The MDH 
activity of the fusion protein was similarly determined at 
340 nm in reaction solutions containing 10 mM malate, 4 mM 
NAD"*", and 0-40 jxm oxaloacetate or 0-50 jxm NADH. 

The citrate synthase (CS) activity of the fusion protein was 
determined using O.I mM acetyl-CoA, 1 mM dithionitro- 
benzoate, and varied concentrations of oxaloacetate (0-1 mM). 
Reactions were monitored at 412 nm, where the reaction 
product of CoA and dithionitrobenzoate shows maximum 
absorption (e = 13 600 M~'cm~'). The conversion of malate 
into citrate in the coupled reaction catalysed by the MDH and 
CS activities of the fusion protein (about 90 U L~^) was 
similarly determined at 412 nm, using 10 mM malate, 4 mM 
NAD"*", 0.1 mM acetyl-CoA, and 0.4 mM dithionitrobenzoate. 
In experiments where the effect of AAT on the coupled reaction 
was examined, 4 mM glutamate and 10-1000 nM AAT were 
added to the reaction solution. 

Kinetic parameter values were determined by nonlinear 
regression analysis based on the least-squares fitting method. 
Numerical solutions of kinetic differential equations were 
obtained using integration procedures of the commercial 
computer program mathematica. 



RESULTS 

Preparation and general properties of the fusion protein 

The fusion protein of the yeast forms of malate dehydrogenase 
(MDH) and citrate synthase (CS) was prepared according to the 
protocol elaborated by Srere and coworkers [12]. The protein 
was purified from 30 g wet E, coli, and only MDH activity was 
measured in the purification steps. Recovered amounts (in 



MDH units) of the fusion protein during purification are given 
in Table 1. 

As previously described [11,12], the fusion protein retains 
malate dehydrogenase (MDH) and citrate synthase (CS) 
activities. These can be assayed separately, or in the coupled 
reaction converting malate into citrate with intermediate 
formation of oxaloacetate (OAA) according to the reaction 
scheme in Fig. 1. The concentration of free OAA in the coupled 
reaction can be estimated, in principle, from the inhibition of 
citrate formation caused by aspartate production from OAA in 
the presence of glutamate and aspartate aminotransferase 
(AAT). The corresponding reaction step is included in Fig. 1 . 

In this investigation, the rate behaviour of the reaction 
system in Fig. 1 was examined under the same conditions as 
those used in the kinetic experiments reported by Lindbladh 
et al. [11], i.e. at 25 "^C in 40 mM phosphate buffer (pH 8.1) 
containing 10 mM malate, 4 mM NAD"*", 0.1 mM acetyl-CoA 
and, when applicable. 4 mM glutamate. 

Malate dehydrogenase activity of the fusion protein 

MDH catalysis of the reversible oxidation of malate by NAD"^ 
is known to be governed by Michaelis-Menten kinetics [16] 
and, at fixed concentrations of the two substrates, should 
conform to the steady-state rate equation 

a(l - ^[OAA][NADH]) 



vx = 



\ + p,[OAA] + p2[NADH] -f 33[OAA][NADH] 



^ = 



1 



[malate][NAD+]A:e, 



(1) 



(2) 



where a stands for the reaction velocity in the absence of 
products, pi, and p3 represent kinetic parameters for 
inhibition of the reaction by the products NADH and OAA. 

Figure 2 shows initial steady-state velocities (vi) 
recorded for the MDH activity of the fusion protein in 
the presence of 10 mM malate, 4 mM NAD"*", and varied 
concentrations of OAA. A fit of Eqn (1) to these data 
for [NADH] = 0 gave a= 87 (± 3) nmoI min"* U"*and 
Pi =0.135 (± 0.0 U) fjLM~'. The parameter P2 was similarly 
determined from initial velocity measurements performed in 
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MOH 

Glutamate 



OAA 
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Fig. 1. Reaction system considered. 



Table 1. Purification of the MDH/CS fusion protein from E. coli cells. 



PurificaUon step 


Protein (mg) 


MDH activity (U) 


Specific activity (U mg ') 


Soni cation 


1600 


19 400 


12 


Nickel -agarose chromatography 


17.2 


1750 


102 


Gel filmition sephacryl S200 


7.3 


1550 


213 
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Fig. 2. MDH activity of the fusion protein. Effect of OAA on the steady- 
state rate {v{) of oxidation of malate (10 mM) by NAD (4 niM) in 40 niM 
phosphate buffer, pH 8.1, catalysed by the fusion protein (92 U L~'). 



the absence of OAA and presence of varied concentrations of 
NADH. This gave P2 = 0.0112 (± 0.0005) |xm"^ The para- 
meter p3 was found to be too small (less than 0.001 m-m"^) to 
have any detectable effect on the reaction kinetics at OAA and 
NADH concentrations below 50 p,M; it was put equal to zero in 
all simulations reported below. 

The denominator term containing q in Eqns (1) and (3) 
reflects the product inhibition contributed by the approach to 
equilibrium conditions. The equilibrium constant foi* the 
oxidation of malate by NAD"^ under the present experimental 
conditions was determined by measurement of the amount of 
NADH formed after catalytic equilibration of a reaction 
mixture containing high known concentrations of malate and 
NAD"^, This gave K^^ = 9.9 x 10"^, corresponding to 
^ = 2.5 X 10^ for 10 mM malate and 4 mM NAD"^. 




500 
[OAA] 



1000 



Fig. 3. CS activity of the fusion protein. Steady-state rate (tp^) of the 
reaction of varied concentrations of OAA with acetyl-CoA (0.1 mM) in 
40 mM phosphate buffer, pH 8.1, catalysed by the fusion protein 
(92 U L~')- 



Mlchaelian parameters for reactions involving OAA 

Figure 3 shows initial steady-state velocities (vj) recorded for 
the CS activity of the fusion protein in the presence of 0.1 mM 
acetyl-CoA and varied concentrations of OAA. As would be 
expected from previous reports [17], the Michaelis-Menten 
equation 

V[OAA] 



V2 - 



(3) 



A:^ + [OAA] 

could be well fitted to the experimental data, giving the 
best-fit estimates V ~ 140 (± 8) nmol-min" "'•U"' and 
11.9 (± 0.8) \LM. 
Parameters defined by the steady-state rate equation 



V3 



fc,,,[OAA]c, 



(4) 



+ [OAA] 

for the action of AAT (active site concentration c^) on OAA in 
the presence of 4 mM gtutamate were similarly determined. 
This gave Jt'ca, = 69 (± 1) s"* and K'^ = 133 (± 6) ^,m. 

Rate behaviour predicted by a free-diffusion mechanism 

As substrates in the reactions now considered are present in 
large excess to the proteins, enzymic species in Fig. I can be 
assumed to be in a steady state over the time-hierarchy where 
the main changes in concentration of nonenzymic reactants 
occur. Provided that reactions are carried out for such a short 
period of time that the high initial concentrations of the 
substrates malate, NAD"*", acetyl-CoA, and glutamate remain 
essentially unchanged, one has to consider only the changes in 
concentration of the reactants OAA and NADH to account for 
the time-course of formation of the products citrate and CoA, 
The time-dependence of the latter four concentration variables 
can be readily expressed if the transfer of OAA occurs by a 
free-diffusion mechanism, and is governed by the differential 
equations 

^^[NADH] 

= -z^i (5) 



dt 



^[OAA] 



■ = Vi — V2 — V2 



d [citrate] d [CoA] 



V2 



(6) 



(7) 



dt dt 

where Vj,, v^, and V2 denote the steady-state rates of the 
reactions catalysed by, respectively, AAT and the MDH and CS 
moieties of the fusion protein. 

Kinetics of the coupled reaction 

Figure 4 shows the time-course of CoA formation in the 
coupled reaction initiated by the addition of catalytic amounts 
of the fusion protein (92 U L~*) to a 40-mM phosphate buffer 
solution (pH 8.1) containing fixed high concentrations of the 
substrates malate (10 mM), NAD"** (4 mM), and acetyl-CoA 
(0.1 mM). The time-dependence expected for concentration 
variables in this reaction (if it proceeds by free diffusion) was 
calculated by numerical solution of Eqns (5-7) for V3 = 0 with 
the assumption that Eqns (1-3) apply, and using the kinetic 
parameter values reported in the previous sections. The results 
are given in Fig. 5. The calculated trajectory for the 
concentration of CoA, multiplied with the absorption coeffi- 
cient (e = 13 600 M~' cm"') for the chromophore used to 
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Fig. 4. Tune-course of the fusion-protein-catalysed coupled reaction. 
Absorbance changes recorded at 4 1 2 nm for the conversion of malate into 
citrate and Co A. Conditions as in Fig. 2, except that reactions were 
performed in the presence of 0.1 mM acetyl-CoA and 0.4 mw dithioni- 
trobenzoate. The dashed curve indicates the absorbance changes expected 
for a free-diffusion mechanism. 
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Fig. 6. Effect of AAT on the steady-state rate of the fusion-protein- 
catalysed coupled reaction. Conditions as in Fig. 4, except that reactions 
were performed in the presence of varied concentrations of AAT and that 
the concentration of fusion protein was 82 U L~'. The curve drawn 
represents reaction rates expected for a free-diffusion mechanism. 



monitor the formation of CoA experimentally, is included in 
Fig. 4 for comparison with the experimentally recorded trace. It 
accounts within experimental precision for the observed time- 
course of CoA formation over the entire reaction time interval 
examined (0-3 min). 



Effect of aspartate aminotransferase 

Figure 6 shows the inhibitory effect of varied concentrations of 
aspartate aminotransferase on the pseudo steady-state rate of 
the coupled reaction catalysed by the fusion protein (82 U L"^) 
under the conditions of the experiment in Fig. 4. The effect 
expected for a free-diffusion mechanism was estimated from 
reaction trajectories calculated by numerical solution of 



Eqns (5-7) with the assumption that Eqns (1-4) apply and 
using the kinetic parameter values reported in the previous 
sections. The results are included in Fig. 6 and show that there 
is an excellent agreement between the calculated steady-state 
rates and those observed experimentally over the entire range of 
AAT concentrations tested (10-1000 nM). 

Furthermore, the transient approach of the inhibited reactions 
towards the pseudo steady-state invariably could be satisfacto- 
rily accounted for in terms of a free-diffusion mechanism. This 
is illustrated in Fig. 7 by example of the results obtained for 
reactions carried out in the presence of 10, 50, and 600 nM 
AAT. Inspection of Fig. 7 shows that the transient lag time 
decreases with increasing concentrations of AAT and eventually 



NADH 




0 12 3 

Time (min) 

Fig. 5. Calculated time-course of reactant concentration changes in the 
fusion-protein-catalysed coupled reaction. Trajectories obtained by 
numerical solution of Eqns (1-7) using parameter values reported in the 
text. 
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Fig. 7. Effect of AAT on the transient phase of the fusion-protein- 
catalysed coupled reaction. Time-course of absorbance changes recorded 
at 412 nm for reactions performed in the presence of 10, 50, and 600 nM 
AAT under the conditions in Fig. 6. Dashed curves show the absorbance 
changes expected for a free-diffusion mechanism, as calculated by 
numerical solution of Eqns (2-8). 
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beconfies immeasurably short. This effect was not detected by 
Lindbladh et ai [11] over the more limited range of AAT 
concentrations tested by them. 



DISCUSSION 

Consistency of the kinetic data with a free-diffusion 
mechanism 

The fusion protein of malate dehydrogenase and citrate 
synthase examined in this investigation is identical with that 
described by Lindbladh et ai [11], except that a His-tag was 
introduced at the N-terminus to facilitate purification of the 
protein [12]. This His-tag modification has been found not to 
affect the kinetic properties of the porcine variant of the fusion 
protein [12] and the kinetic results now obtained with the yeast 
variant do not differ in any essential respect from those 
previously reported. The main new information provided here 
does not concern the actual rate behaviour of the fusion protein, 
but the mechanistic interpretation of the rate behaviour. 

Lindbladh et al. [11], in their analysis of the kinetics of the 
coupled reaction catalysed by the fusion protein, assumed that 
the expected time-course for a free-diffusion reaction should 
agree with that observed for an adequately composed mixture 
of the free enzymes; the differences in kinetic properties of the 
free enzymes and the corresponding moieties of the fusion 
protein were considered to be small enough to justify such an 
approach. In this investigation, we have taken the more direct 
approach of basing all predictions of the rate behaviour of the 
fusion protein on the kinetic properties of the fusion protein 
itself. This means that the expected kinetics of the reaction 
system in Fig. 1 can be unambiguously established by solution 
of Eqns (5-7) in case a free-diffusion mechanism applies. 

The results in Figs 4 and 5 therefore provide the indisputable 
inference that the first 3 min of the observed time-course of the 
fusion-protein-catalysed malate to citrate conversion can be 
excellently accounted for with the assumption that the reaction 
conforms to the scheme in Fig. 1 with OAA transfer by free 
diffusion. This observation has direct bearing on the mechan- 
istic problem now considered; conclusions drawn by Lindbladh 
et aL [11] were based on the behaviour of the fusion protein in 
the reaction phase reflecting the transient approach to steady- 
state conditions, and this reaction phase is well covered by the 
data in Fig. 4. 

Arguments previously presented in support of channelling 

Lindbladh et al. [11] argued that transient lag times for the 
coupled reaction catalysed by the fusion protein are indicative 
of OAA channelling because they are much shorter than 
expected for a free-diffusion mechanism. The results in Fig. 4 
invalidate this argument by showing that the experimental trace 
during the first 3 min of the reaction differs insignificantly from 
that calculated for a free-diffusion mechanism. This means that 
the observed transient time, as well as the observed steady-state 
rate of the reaction, agrees with that predicted by the kinetic 
properties of the MDH and CS moieties of the fusion protein. 
As the transient time and steady-state rate both reflect the 
steady-state level of OAA that is reached after the transient 
reaction phase, data in Fig. 4 provide clear evidence that OAA 
attains a steady-state concentration agreeing with that expected 
for a free-diffusion mechanism. The data lend no support to the 
idea that the fusion protein acts upon a catalytically effective 
concentration of the intermediate OAA that is higher than the 



concentration of OAA in bulk solution or otherwise leading to a 
more rapid transfer of OAA than that obtained by free diffusion. 

In view of these results, the observation [II] that the transient 
time for the coupled reaction is longer with the free enzymes 
than with the fusion protein looses relevance with regard to the 
mechanism of action of the fusion protein and merely indicates 
that a higher steady-state level of OAA is reached with the free 
enzymes. Such a difference in behaviour of the two systems 
seems to be fully consistent with the reported differences in 
ATmCOAA) values for the MDH and CS activities of the free and 
fused enzymes [11]. The net balance of these value 
differences is in the direction that would be expected to favour 
a lower steady-state level of OAA in the fusion protein. 

Similar considerations apply for the effect of AAT on the 
fusion-protein-catalysed coupled reaction. Data in Figs 6,7 
establish that the experimentally determined effects of AAT on 
the steady-state and presteady-state time-course of the reaction 
are in satisfactory agreement with those expected for a free- 
diffusion mechanism according to the rate parameter estimates 
obtained for the enzymic activities involved. Consequently, 
there is no longer any reason to believe that the fusion protein 
sequesters the intermediate OAA and prevents it from being 
fully accessible to AAT. The observation that the coupled 
reaction is more strongly inhibited by AAT when catalysed by 
the free enzymes is rendered irrelevant by the results in Figs 6,7 
and has to be explained in terms of the existing differences in 
kinetic parameter values for the free and fused enzymes. 

General conclusions 

Summing up the above discussion, the present results invalidate 
the main kinetic arguments that led Lindbladh et aL [11] to 
propose that there is a channelled transfer of the intermediate 
OAA in the coupled reaction catalysed by the fusion protein. It 
may now be concluded that the rate behaviour of the fusion 
protein of yeast MDH and CS is fully consistent with that 
expected for a free-diffusion mechanism of OAA transfer. This 
finding provides some additional inferences of general interest. 

First, ambiguities obviously arise when the kinetics of 
coupled reactions catalysed by fusion proteins are interpreted 
using the rate behaviour of the free enzymes as a predictive 
model of the expected free-diffusion behaviour of the fused 
enzymes. All claims so far made for the kinetic detection of 
metabolite channelling in fusion proteins would seem to have 
been based on such an interpretational approach. There is 
reason to consider these claims as tentative until they have been 
confirmed by less ambiguous methods such as the one 
described in this investigation. 

Second, bringing two sequential enzymes together in a fusion 
protein (or, by inference, in a bienzyme complex or through 
adjacent attachment to cell structures) is obviously not 
sufficient to cause any kinetically significant metabolite 
channelling through proximity effects. This is consistent with 
the conclusion drawn by Elcock & McCammon [18] from 
Brownian dynamics simulation studies of a tentative model 
structure of the fusion protein of MDH and CS, which indicated 
that less than 1% of the OAA molecules produced at the MDH 
site would be directly (i.e. without prior release to solution) 
transferred to the CS site through Brownian motion alone. 

However, Elcock & McCammon [18] also found that a 
positive electrostatic surface exists between the two enzymic 
sites, and that the inclusion of electrostatic interactions in the 
Brownian dynamics simulations increased the estimated direct 
transfer efficiency to around 45%. This lead them to propose 
that electrostatic surface diffusion accounts for a highly 
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efficient channelling of OAA in the fusion protein, and 
subsequent theoretical work [19] indicated that this proposal 
could be reconciled with the kinetic effects that Lindbladh et al. 
[II] attributed to channelling. The present results call for a 
revaluation of the calculations performed by Elcock and 
collaborators [18,19] and lend no support to their proposal. 
The data in Figs 4-7 establish that if surface diffusion con- 
tributes to the transfer of OAA. then these contributions are too 
small to bring the MDH-produced OAA out of apparent 
equilibrium with bulk solution and to cause any kinetically 
significant deviations from the rate behaviour predicted by a 
free-diffusion mechanism. 

Finally, it should be mentioned that Eqns (5-7) were derived 
with the tacit assumption that the MDH and CS moieties of the 
fusion protein operate independently of each other, such that 
kinetic parameters for the MDH activity remain unaffected by 
catalytic events at the CS site (and vice versa). The validity of 
this noncooperativity assumption is corroborated by the 
agreement between calculated and observed data in Figs 4-7. 
The only detectable kinetic consequence of the fusion of the 
two enzymes therefore is the previously reported modification 
of the kinetic parameter values for the MDH and CS activities 
from those of the free enzymes to those of the fusion protein 
[II]. These parameter modifications (and the underlying rate 
constant changes) could be of biological interest as an illus- 
tration of effects that might be caused by complex formation 
between consecutive enzymes in metabolic pathways. For 
example, and by analogy to what seems to be the case for the 
fusion protein of MDH and CS, changes in the kinetic 
parameter values could result in lower steady-state levels of 
accumulating intermediates and hence in shorter transient times 
for the response of the pathway to changes in the environmental 
conditions. Such parameter-related effects have nothing to do 
with metabolite channelling, however, but can be adequately 
described and understood in terms of metabolite transfer by free 
diffusion. 
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